The role of weak charging in metastable colloidal clusters 
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We study metastable clusters in a colloidal system with competing interactions. A short-ranged 
polymer-induced attraction drives clustering, while a weak, long-ranged electrostatic repulsion pre- 
vents extensive aggregation. We compare experimental yields of cluster structures expected from 
theory, which assumes simple addition of the competing isotropic interactions. For clusters of size 
4 < m < 6, the yield is significantly less than that expected. We attribute this to an anisotropic 
self-organized surface charge distribution linked to the cluster symmetry: non-additivity of electro- 
static repulsion and polymer-induced attraction. 7-membered clusters have a clear optimal yield of 
the expected pentagonal bipyramid structure as a function of strength of the attractive interaction. 
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Clusters are intrinsically distinct from bulk matter, 
with the restricted degrees of freedom making their en- 
ergy landscapes tractable, providing ideal model systems 
for studying kinetic pathways [J . Although atomic clus- 
ters have been extensively studied, direct visualization is 
challenging and often limited to states of low temperature 
|2 . Colloidal dispersions provide a new medium by which 
to study clusters, and their well-defined thermodynamic 
temperature allows analogy with atomic, molecular, and 
protein clusters. Furthermore, the ability to track parti- 
cle coordinates in 3D yields detailed static and dynamic 
information usually only accessible to computer simula- 
tion [3]. 

Colloidal dispersions feature novel cluster phases [4, 
[3 m [7] , which open the possibility to synthesize colloidal 
clusters and 'molecules' [8j |9l [TO] HTJ [12^ presenting excit- 
ing possibilities for new materials [13] . To fully realize the 
potential of colloidal clusters and 'molecules', it is nec- 
essary to understand their behavior, and in particular to 
investigate ways to optimize yields of desired structures 
[9]. Perhaps the simplest cluster- forming system is built 
around a spherically symmetric attraction which may be 
realized by adding non-adsorbing polymer to a colloidal 
suspension. The polymer osmotic pressure then leads to 
an effective attraction between the colloids known as the 
Asakura-Oosawa (AO) potential, whose strength is ap- 
proximately proportional to the polymer concentration 
Cp [Fig. [Ija)]. 

While colloidal systems with purely attractive interac- 
tions may aggregate, stabilization can be achieved with 
longer-ranged electrostatic repulsions [5 , which are of- 
ten treated with a Yukawa (YUK) potential [14 . Here 
we combine the competing short-ranged attractions and 
weak long-ranged repulsions in colloidal systems with 
the well-understood behavior of clusters interacting via 
short-ranged attractions [T5^. We neglect possible many- 
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Figure 1: (color online) (a) Various potentials: Morse (black 
line) denotes a Morse potential for a well depth £m = 4.6 and 
pQ = 33.1. AO (red dashed line) denotes Asakura-Oosawa 
for Cp= 6.02 X 10"^, YUK (light cyan line) denotes a Yukawa 
interaction for a contact potential of k^T and inverse Debye 
screening length = 0.5. AO+YUK (blue dotted line) sums 
AO and YUK illustrates a treatment of the electrostatic re- 
pulsions as a small perturbation (see text), (b) Determining 
the repulsive interactions between the clusters. Experimental 
g(r) for Cp = 6.02 x 10~^ (circles) and MC simulation (solid 
line) treating clusters as single particles, each with a charge 
Z = 61. (c) Confocal micrograph with Cp = 1.43 x 10""^. Bar 
is 10 /im. 



body interactions [16 . As Fig. [TJa) shows, for short in- 
teraction ranges the AO potential is similar to the Morse 
potential, whose range is set by a parameter po- Now 
the structures of the global energy minima for clusters of 
Morse particles are known [15], and for small clusters of 
size m < 8, the topology of this minimum structure is 
not sensitive to the range of the interaction, nor to weak, 
long-ranged Yukawa repulsions [18 . With the exception 
of m = 6 (see below), the global minimum of small clus- 
ters in our colloidal system should therefore be identical 
to those tabulated for the Morse interaction [15 . Us- 
ing a relatively dilute colloidal dispersion, we consider 
each cluster as a separate subsystem, and investigate the 
cluster structure at the single particle level. Using the 
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topological cluster classification algorithm [17], we com- 
pare the cluster structures found with those correspond- 
ing to Morse global energy minima, as a function of poly- 
mer concentration (strength of the attractive interaction) 
[151 llfii- As we shall see, even in our simple model sys- 
tem, the behavior is not as straightforward as one might 
expect. We believe this is due to a breakdown of the 
Yukawa description of the electrostatic interactions. 

We used poly(methyl methacrylate) (PMMA) colloids 
sterically stabilized with polyhydroxyl steric acid. The 
colloids were labeled with fluorescent rhodamine dye and 
had a diameter of <j = 2.0 /im and a polydispersity of 
4% as determined with static light scattering. The van 
der Waals interactions are neglected by closely matching 
the refractive index between the solvent and particles, 
sedimentation is avoided by matching the density. The 
colloid volume fraction was (j) ^ 0.02. We use a solvent 
composition of approximately 0.373 cis decalin, 0.093 cy- 
clohexyl bromide and 0.533 tetrachloroethylene by mass, 
in which the dielectric constant ~ 2.2. This is a good 
solvent for the polymer, polystyrene (PS). We used a 
molecular weis; ht of 3.0 X 10^, where our measurements 
of the interactions agreed well with the Asakura-Oosawa 
(AO) model with some swelling (35%) of the polymer 
[19], leading to a polymer-colloid size ratio q ^ 0.22 which 
corresponds to a Morse potential with po = 33.1 [2^. We 
studied the system using confocal microscopy at the sin- 
gle particle level, using a Leica SP5, with which we track 
the particle coordinates to an accuracy of around 100 
nm. We impose a bond length of 1.4cr, this slightly larger 
value than the range of the attractive interaction reflects 
particle tracking errors and polydispersity. Our results 
are robust to reasonable changes in the bond length. We 
sampled typically 10^ coordinates per state point. No 
change in the cluster populations was observed on the 
timescale of the experiments (1 day). 

We find the order of magnitude of the colloid charge 
in cluster fluid by treating the clusters as individual par- 
ticles and neglecting their different sizes [21]. The size 
distribution is shown in Fig. |2|a) for a polymer weight 
fraction Cp = 6.02 x 10~^ the mean < m >= 3.3 and 
standard deviation is 4.4. In Fig. [ijb), we fit the ra- 
dial distribution function g{r) with Monte Carlo sim- 
ulation according to a Yukawa potential [14 , yielding 
Z = 61 ± 20 charges per cluster, and a Debye screening 
length of K,~^ = 2.0 ± 0.4cr. This suggests around 10 — 20 
charges per colloid, which corresponds to a Yukawa po- 
tential at contact of Psyvk = ^^/B/[(l+^cr/2)^cr] ~ 1 — 3 
where /3 = 1//cbT (/cb: Boltzmann's constant) Ib ^ 25 
nm is the Bjerrum length. Although this charge is very 
low, similar values have been measured for PMMA in 
apolar solvents [22] . A confocal microscopy image of the 
system is shown in Fig. [ijc), we see the clusters are well 
isolated, consistent with our assumption that each clus- 
ter may be treated as a separate system. Increasing Cp 
promotes clustering. 

We now proceed to consider the behavior of the indi- 
vidual clusters (as a separate system). Figure [2|b) shows 
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Figure 2: (color online) (a) Cluster size distribution for Cp= 
6.02 X lO"'^. (b) 'Structural yield' Nc/Nm as a function of Cp. 
Only in the case of m = 3 does the yield of the global energy 
minimum (3 A) approach 100%. Vertical dashed line denotes 
polymer overlap. ISsao is the contact potential for the AO 
interaction. 



the 'structural yield' Nc/Nm as a function of the attrac- 
tive interaction Cp. Nc is the number of clusters of a 
given structure, N^ is the total number of clusters con- 
taining m colloids. For the smaller cluster sizes a greater 
portion are found in the Morse global minimum, with 
the exception of the octahedal 6A cluster which is much 
less popular than either the triangular bipyramid 5A or 
the pentagonal bipyramid 7A. Considering the topology 
of the 3A, 4A, 5A, and 6A clusters we see that the dis- 
tance between all the nearest neighbors will be the same 
due to the short range nature of the AO potential and 
next nearest neighbors do not contribute, if we ignore the 
weak electrostatic repulsion. The minimum energy for 
these clusters is thus approximately proportional to the 
number of nearest neighbors or contacts. However, for 
m = 6, global minimum for the Dzugutov potential [23] 
has a C2v point group symmetry and the same number 
of contacts as 6A. We therefore seek the this 6Z cluster 
using the TCC. 

Neglecting the weak charge, both 6A and 6Z clusters 
have the same ground state potential energy but there are 
other contributions to the cluster's free energy: configu- 
rational entropy and electrostatic energy. Indistinguish- 
ablility arguments suggest that 6Z should be 12 times 
as popular as 6A. Moreover, 6Z is further favored by its 
larger radius of gyration, which lowers electrostatic en- 
ergy [24 . As we shall see below, the effects even of these 
weak electrostatic interactions can be rather complex. 
However, the very large population difference in this sys- 
tem is suggestive of a further mechanism at play: the 
kinetic pathway and the associated energy barrier. Addi- 
tion of particles leads to the sequence 3A^4A^5A^6Z 
while the formation of 6A requires a bond to be broken 
and insertion of a particle to form the 4-membered square 
ring. 

We expect an increase of polymer concentration (low- 
ering of the effective temperature) to promote clusters 
which maximize the number of bonds, i.e. 3-7A, or 6Z. 
This we see in Fig. |2|b) for m = 3, in the near 100% 
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Figure 3: (color online) (a) Distribution of bond angles Q 
for m = 3 clusters. As Cp (listed xlO~^) is increased, the 
population of 3 A triangles increases, (b) Maximum center 
separation, m = 4. Red dashed line denotes 4A tetrahedron, 
blue (solid) line 4D diamond and orange (dotted) line ('4X') 



higher energy (fewer bonds). 
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yield of 3 A triangles at high Cp. In Fig. [s] (a) we show 
the distribution of bond angles ^ in m = 3 at different 
polymer concentrations. In addition to a higher 3 A trian- 
gle population at higher Cp, we see a peak at large bond 
angle, showing that an elongated structure is preferred, 
and that the system seems to sit either in a triangle or the 
linear structure. We attribute this two-state behavior to 
the fact that the colloids carry a (weak) charge, promot- 
ing elongated states [24 . In the case of more complex 
clusters, such as 7A, an initial increase in yield with Cp 
is followed by a decrease, presumably due to frustration, 
at high interaction strength. 

However, intermediate cluster sizes, 4 < m < 6 show 
a surprising behavior, with little sensitivity to polymer 
concentration. Motivated by our analysis of m = 3 clus- 
ters, we consider structures in addition to the ground 
state 4A tetrahedron in the m = 4 population. Here 
we consider the largest separation of colloid centers d. 
Now d <j for 4A tetrahedra, while the maximum value 
is d 3cr for a line of 4 particles. Figure |3jb) shows 
two dominant structures in the m = 4 population: the 
4A tetrahedra and a diamond-shaped structure we term 
4D, which is distinct from a square structure for which 
d = \/2a. Unlike more complex (e.g. 7A) clusters, for 
m = 4, there is no geometric frustration in the forma- 
tion of tetrahedra. In fact Brownian dynamics simula- 
tions of four particles interacting via the Morse potential 
[25 show a clear peak in the 4D population, and a ten- 
dency to form tetrahedra at strong interaction strength 
[Fig. [4]^ a)]. Furthermore, using an isotropic Yukawa in- 
teraction with parameters similar to those found in our 
experiments has little effect: we see a high yield of 4 A at 
large /3£m^ consistent with [18 . 

Why then do the experiments have such a high 4D 
population? One argument might be that the AO model 
is starting to break down as we exceed overlap, which oc- 
curs at Cp ^ 8.5 X 10 ~^ as indicated in Fig. 2] However 
for this size ratio we expect only a slight effect in the 
AO interactions at overlap [26 . Instead, we postulate 
that 4D is stabilized by the weak charging present in this 
system, which unlike the spherically symmetric distribu- 
tion often assumed, may be anisotropically distributed 
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Figure 4: (color online) (a) Structural yields Nc/Nm as 
a function of well depth /3£m calculated from Brownian dy- 
namics simulations for the Morse potential for po = 33.1 
for m — 4. Diamonds denote 4D and triangles 4A. Thick 
lines with unfilled symbols denote Yukawa contact potential 
Psyvk = 0, light blue (dashed) Psyvk = 1 and thin crim- 
son lines and filled symbols pevvK = 3 and inverse Debye 
screening length = 0.5. (b) Conductivity of the solvent, 
and colloidal suspension at various polymer concentrations 
xlO"^ 



on colloid surfaces, leading to anisotropic electrostatic 
interactions [7 . 

While electrostatic charging mechanisms in these non- 
aqueous systems are not yet fully understood, it is reason- 
able to suppose that charging may be related to the free 
colloid surface, which may be reduced upon clustering. 
To this end, we measured the conductivity of the sys- 
tem as a function of Cp with dielectric spectroscopy [Fig. 
|4jb)]. As we do not add salt, we assume the system is in 
the counter-ion dominated regime, and we note a large in- 
crease in conductivity relative to the pure solvent. What 
is evident from Fig. |4jb) is that the introduction of the 
colloids massively increases the conductivity relative to 
the pure solvent, and that adding polymer then decreases 
the conductivity. We associate this reduction in ionic 
strength (charge per colloid) with clustering, and argue 
that the charge is indeed associated with the free surface 
of the colloids, which is reduced upon clustering. While 
precise determination of the ionic strength from conduc- 
tivity measurements requires detailed knowledge of the 
ionic species present, assuming Walden's rule and a typ- 
ical limiting molar conductance for ions in these systems 
of 4 cm^mol~^fT4^ we arrive at an ionic strength around 
10~^^ M, consistent with the Debye length K:~^ ^ 2cf 
found with our MC simulations [Fig. [ijb)]. 

4D and 4A clusters are formed by an attachment of one 
particle to 3A. For 3A^4A there are three direct contact 
points to be formed, whereas for 3A^4D there are only 
two. Thus 4D may be more easily formed than 4A due to 
a reduced electrostatic barrier and smaller loss of trans- 
lational entropy of counter ions. We note that the latter 
factor is not considered at all in a model with a Yukawa 
potential. Furthermore, we stress that a Yukawa poten- 
tial should break down at short distances and thus fail 
in describing aggregation of particles. For 4D (i ~ >/3cr 
whereas for 4 A d ^ a. So 4D is more favored than 4 A 
electrostatically. We note that although there is no struc- 
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tural frustration in the 4D^4A transition, there is only 
one path, if no bonds are broken; one or both of the 
two end particles must 'roll' around the central two. We 
note that the 4D-4A transformation, always leads to a 
decrease in which may be prohibited by electrostatic 
repulsions. Thus, once 4D is formed, the transforma- 
tion from 4D to 4 A is suppressed: kinetic arrest. That 
4D formation is promoted and 4D^4A transformation is 
suppressed due may lead to the high fraction of 4D. It is 
possible that the charging mechanism leads to a scenario 
in which the equilibrium populations are a distribution 
dominated by both 4D and 4A as found. However, the 
lack of response of the system to increasing the strength 
of attraction by a factor of three leads us to the con- 
clusion that some form of kinetic arrest is more likely. 
On the other hand, our results suggest that it is not so 
difficult to to make a 3 A triangle from a linear cluster. 
Although the reason for the difference is not immediately 
clear, we note that the linear^3A transition is not geo- 
metrically constrained to follow a single path. 

Here we consider some microscopic aspects of this 
problem. The low surface charge density (small Z) im- 
plies that the discrete nature of charges may be impor- 
tant. Furthermore, the recombination rates of ions might 
be very slow: even for small ions, diffusion-controlled re- 
action rate constants are 10^-10^^ dm^M^^s""*^, of order 
seconds or tens of seconds here. For the larger ions or 



complexes relevant to these apolar solvents [T4[ |22] the 
rate may be lower still. That the charges are discrete in 
both space and time may play an important role in the 
above-described frustration. Further work is necessary 
for elucidating the exact mechanism responsible for the 
observed cluster type distribution. 

Our study suggests that the one-component descrip- 
tion is breaking down, even for these simple colloidal 
systems. In the case of m = 6, two structures with simi- 
lar energies compete, however, the 6Z is very much more 
popular than the 6A octahedron. Although equilibrium 
arguments favor this trend we note that 6Z is structurally 
similar to 5 A, and thus may likely be kinetically favored. 
The yield of the fivefold symmetric 7A shows signs of frus- 
tration at relatively deep quenches, reminiscent of patchy 
particles [9 . Our results form the beginnings of an un- 
derstanding of self-assembly of colloidal molecules as a 
function of attraction strength; we hope they will stimu- 
late the development of more powerful models. 
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